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Characterization of the high-pressure superconductivity in the Pnma phase of calcium 
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The thermodynamic parameters of the superconducting state in calcium under the pressure at 161 
GPa have been calculated within the framework of the Eliashberg approach. It has been shown that 
the value of the Coulomb pseudopotential is high = 0.24) and the critical temperature {Tc ~ 25 
K) should be determined from the modified Allen-Dynes formula. In addition, it has been found that 
the basic dimensionless ratios of the thermodynamic parameters significantly diverge from the BCS 
predictions, and take the following values: (i) The zero temperature energy gap to the critical tem- 
perature {Ri = 2A (0) /fcsTc) is equal to 4.01. (n) The ratio R2 = (C^ (Tc) - {Tc)) /C^ (Tc) 
equals 2.17, where and denote the specific heats for the superconducting and normal state, 
respectively, (iii) The quantity R3 = TcC^ {Tc) /He (0) = 0.158, where He indicates the thermo- 
dynamic critical field. Finally, it has been proven that the electron effective mass is large and takes 
the maximum of 2.32me at Tc. 
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I. INTRODUCTION 

Under the influence of the high pressure, the elemental 
calcium undergoes a series of structural phase transitions. 
In particular, one can distinguish seven phases in the 
range of pressure {p) from to 241 GPa [H, 0, @ (please 
see Figure [T] (A) for the details). 

The two first phases, namely Ca-I and Ca-II, have been 
classified as a fee and bee structures, respectively [H, 
0. The third phase (Ca-III) has been primarily linked 
with the sc structure, however the recent reports suggest 
other assignments. On the basis of the theoretical stud- 
ies, Teweldeberhan et al. proposed the Cmmm structure 
Nakamoto et al. also vote in favor of the Cmmm 
structure @. On the other hand, Mao et al. have pre- 
dicted the transition from the sc-like structure to the 
monoclinic phase at 30 K and p ~ 40 GPa Q. It needs 
to be underlined that the stability of the structure sc in 
the area of the existence of the phase Ca-III is being con- 
firmed by the results achieved by Errea et al. and Yao et 
al, at least for the temperature of 300 K 0, @|. 

The existence of the phases Ca-IV and Ca-V has been 
experimentally examined in the papers Q and [H. Fu- 
jihisa et al. have proposed for them the following as- 
signment: the structure Ca-IV should be characterized 
by P4i2i2 and Ca-V by Cmca space groups, respectively 
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FIG. 1: (A) The sequence of the structural phase transitions 
in calcium determined on the basis of the experimental data. 
(B) The dependence of the critical temperature on the pres- 
sure: stars - Okada et al. [l^li squares - Yabuuchi et al. [l5| . 
circles - Sakata et al. 



In the year 2010, Nakamoto et al. discovered the new 
Ca-VI phase with the Pnma structure [l3|- Further, in 
the year 2011, Sakata et al. have reported the existence 
of the host-guest phase Ca-VII We can notice that 
the high-pressure phase of the host-guest character had 
been previously predicted by Arapan et al. and then by 
Ishikawa et al. [ul, [Hi. 

The first mention of the existence of the pressure- 
induced superconducting state in calcium was provided 
by Dunn and Bundy in 1981 [l^. Fifteen years later, 
Okada et al. determined the dependence of the critical 
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temperature (Tb) on the pressure up to the 150 GPa [Tj| 
(please see Figure [1] (B) for details). In the year 2006, 
Yabuuchi et al. have repeated the experimental studies 
of Okada [l5j. It has been found that the values of the 
critical temperature increase much faster together with 
the increase of p in comparison to the results achieved by 
Okada. The last notable experimental results have been 
obtained by Sakata et al. [J] . On the basis of Figure [1] 
(B), it can be easily noticed that for p ~ 216 GPa, the 
critical temperature takes the value equal to 29 K (the 
highest observed Tc among all elements). However, this 
result has been challenged by Andersson [l6| . 

In the presented paper, we have determined all rele- 
vant thermodynamic parameters of the superconducting 
state that is induced in calcium under the pressure at 161 
GPa. We draw the readers' attention to the fact that the 
pressure of 161 GPa represents the highest value of p con- 
sidered by Yabuuchi et al. [l5|. Additionally, the high 
value of the critical temperature at p = 161 GPa. which 
is equal to ^ 25 K, has been recently confirmed by the 
results obtained by Sakata et al. 0. 

For the purpose of this paper, we have assumed that 
the phase Ca-VI is being characterized by the Pnma 
crystal structure. To support this assumption we quote 
the results presented in: [l^l, [l3l and fl8| . 



II. THE ELIASHBERG EQUATIONS 

On the imaginary axis {i = 1), the order parameter 
(A„ = A (iujn)) and the wave function renormalization 
factor {Zn = Z (iu>n)) can be calculated by using the 
Eliashberg equations [l^l : 
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In equations ([T|) and ([2]) the symbol a;„ = ^ (2?i — 1) de- 
notes the 71-th Matsubara frequency, where /3 = l/ksT, 
and kB is the Boltzmann constant. 

The complicated form of the electron-phonon pairing 
kernel is represented by the expression: 



K (Wn - Wm) 
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{uJn - Wm) + 
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where fimax is the maximum phonon frequency (ilmax = 
71.37 mcV); and a^F (fi) indicates the Eliashberg func- 
tion, which models the shape of the electron-phonon in- 
teraction in a detailed way. In the presented paper, the 
form of the a^F (51) fimction has been taken from Yin et 
al. [13. 

The depairing interaction between electrons is de- 
scribed with the use of the Coulomb pscudopotential 



(/i*). The symbol 9 denotes the Heaviside function and 
UJc is the phonon cut-off frequency: lOc = Silmax- 

In the presented paper the Eliashberg equations have 
been solved for 2201 Matsubara frequencies (M = 1100). 
In this case, the obtained solutions are stable for the tem- 
peratures greater than or equal to Tq = 5 K. A detailed 
discussion of the numerical method has been presented 
in [20]- ;24|. 



III. THE COULOMB PSEUDOPOTENTIAL 

The physical value of the Coulomb pscudopoten- 
tial (fi^) can be defined by using the condition: 



[A,, 



= 1\t=Tc 



0, where the critical temperature is equal 



to the experimental value {Tc = 25 K) [15| 
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FIG. 2: (A) The dependence of the order parameter on ojm 
for selected values of the Coulomb pseudopotential {T = Tc). 
(B) The maximum value of the order parameter as a function 
of the Coulomb pseudopotential. 

In Figure [2] (A) we have presented the dependence of 
the order parameter on Um for selected values of 
One can notice that together with the increase of the 
Coulomb pseudopotential, the largest value of the order 
parameter {Ajn=i) decreases. Additionally, in Figure [2] 
(B) we have outlined the complete form of the function 
A,„-i (//*). On the basis of the obtained results, we have 
found that the physical value of the Coulomb pseudopo- 
tential is equal to 0.24. 

The above result means that the depairing electron 
correlations in calcium are relatively strong (the classical 
low-tcmpcraturc superconductors fi^ is about 0.1 [25j). 
It can be noted that similar non-standard value ofu^ has 
been obtained for lithium and CaLi2 ^22,], [l^, [29j . 
[sot ■ For example, the properties of the superconducting 
state in the fee phase of lithium for the pressure values 
22.3 GPa (Tc = 7.27 K) and 29.7 GPa {Tc = 13.93 K) 
have been specified in the paper (29j . It has been shown 
that the physical value of the Coulomb pscudopotential 
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increases with p from 0.22 to 0.36. In the case of CaLi2, 
the parameter is equal to 0.23 {p = 45 GPa and Tc = 
12.9 K) [13]. 

The high values of /x^ in Ca, Li, and CaLi2 are 
difficult to explain in the framework of the classi- 
cal Morel- Anderson (MA) model [3l|, where: = 
/i [1 -)- /xln (we/wpTi)] ^. The symbol /i is defined by: 
fi = p{0)Uc, where p{0) is the electronic density of 
states at the Fermi level and Uc is the Coulomb poten- 
tial; u!e and Ljph denote the characteristic electron and 
phonon frequency, respectively. Since Ue ^ Wp/i, the MA 
pseudopotential is of the order 0.1 and /x^ ^ fi. 

It can be noted that the MA model corresponds to 
treating the irreducible vertex to the first order in Uc- 
Recently, Bauer, Han, and Gunnarsson have extended 
the MA theory to the second order in Uc- The main 
result is that the retardation effects lead to the reduction 
fi — >■ ^J. also in the higher order calculation, but not as 
efficiently as in the first order (s^ . The model presented 
in the paper (s^ is probably the most advanced attempt 
to explain the high values of in Ca, Li, and CaLi2. 

In particular, Bauer et al. have given the following 
expression for the physical value of the Coulomb pseu- 
dopotential: 



fJ-c = 
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where a — 1.38 and a ~ 0.10. 

On the basis of the equation ([4]) one can easily estimate 
the value of Uc for the real materials. In the paper, we as- 
sume the following: oje = W {W is the half-band width). 



= UJln (win = CXp 



f/„"-df7^1n(n)J), and 

p{0) = 1/2 (the constant DOS). In the literature, the 
values of all the important parameters are provided only 
for CaLi2 {p = 45 GPa). In particular: p % = 0.23, 
W = 1991 meV, and wi„ = 17.02 meV [Ig, lH. The 
result is the following: Uc — 2803 meV. 

Next, we address an important issue, namely, how big 
is the error bar of the calculated physical value of the 
Coulomb pseudopotential (A/i^). First of all, we can 
notice that in the framework of the presented analysis, 
the value of p^ depends on the shape of the Eliashberg 
function and the accuracy of the experimental value of 
Tc- For calcium the appropriate Eliashberg function, 
taken from [l3], has been calculated by using the linear- 
response method (full-potential LMTART code [3i|). In 
that paper, the Eliashberg function error bar has been 
omitted ([A/i^]^2j7(n) ~ other hand, the value 

of the critical temperature has been measured with the 
accuracy about ±1 K [l^. On the basis of these facts, 
we have obtained: [A/iJ,]^, = ±0.02. 



IV. THE CRITICAL TEMPERATURE 
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FIG. 3: The dependence of the critical temperature on the 
Coulomb pseudopotential. The filled circles mean the re- 
sults obtained by using the Eliashberg equations; the arrow 
indicates the experimental value of the critical temperature 
{fiQ — 0.24). The solid line denotes the results obtained with 
the help of the modified Allen-Dynes formula. Finally, the 
dashed and dotted lines have been generated based on the 
classical Allen-Dynes formula and the McMillan expression, 
respectively. 



the basis of Eliashberg equations. However, in the case of 
the data interpretation, it is far more convenient to use 
the simple formula that explicitly reproduces the results 
of the advanced numerical calculations. In the branch 
literature there arc known two basic formulas that serve 
for the determination of the critical temperature's value. 
The first one has been introduced by McMillan [3^; the 
second one is the Allen-Dynes expression [s^. Unfortu- 
nately, in the case of calcium, both formulas consider- 
ably underestimate the critical temperature. According 
to the above, the Allen-Dynes formula has been modified 
in such a way, that it allows us to reproduce the numeri- 
cal results correctly. Particularly, in order to achieve the 
proper values of the fitting parameters, the dependence of 
the critical temperature on the Coulomb pseudopotential 
has been analyzed on the level of the Eliashberg equations 
(only a^F (il) has been considered as the physical input 
parameter). Next, the least-squares method was applied. 
The obtained result is presented below: 



ksTc = fi [p*) h (m*) 
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where the functions /i (//*) and /2 (/i*) are expressed by 
the formulas: 



In the framework of the presented formalism, the exact 
value of the critical temperature should be obtained on 
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The parameters, that depend on the Ehashberg func- 
tion, can be determined on the basis of the expressions: 
UJ2 = J /o""^^ dfla^F (n) and A = 2 /g""""" dn^Lmi, 
For calcium under the pressure at 161 GPa, we have re- 
spectively: y^LJ2 = 34.36 meV and A = 1.27. The fitting 
functions Ai (/i*) and A2 (/i*) are presented in the follow- 
ing forms: 



and 



Ai EE 0.145(1 + 115.862/i*), 



A2 (a**) = 5.185(1 - 2.247Ai*) 
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In Figure [3] we have presented the numerical solutions 
obtained with the use of the Eliashberg equations and the 
modified Allen-Dynes formula. Additionally, for the com- 
parison purposes, we have depicted the results based on 
the classical formulas derived by Allen-Dynes and McMil- 
lan. On the basis of Figure [31 one can observe, that the 
expression (4) perfectly reproduces the exact Eliashberg 
numerical predictions. 

The constants in the equation (4) deviate notably 
from the original parameterization. This situation is 
connected with the fact that the analysis based on 
the real-axis Eliashberg equations suggests only the 
semiphenomenological form of the T^-formula: ksTc = 



■ exp 



-b(l+A) 



(see the detailed discussion in [2£ 



X-fi* {l+cX) 

p. 1051-1052). In the case of calcium under the pressure 
at 161 GPa, the values of the Allen-Dynes parameters 
are inappropriate. Thus, the constants (a ^ c) should 
be fit to the data taken from the exact solutions of the 
Eliashberg equations. We can notice that the change of 
the coefficient in the expression (4) under uji^ from 1.2 
to 1.45 slightly lowers the phonon frequency; the two re- 
maining parameters (1.03 and 0.06), in comparison to 
the classical parameterization, increase the value of the 
effective electron-phonon coupling constant. 

Moreover, the parameterization of the strong-coupling 
correction function (/i (fj,*)) and the shape correction 
function (/2 (fJ-*)) also deviates from the original form. 
The achieved result indicates that for the high-pressure 
superconducting state in calcium the shape function has 
greater significance than in the classical superconductors. 

The value of the critical temperature for p = 160 GPa 
has been also calculated in the paper (T^]. By using the 
Allen-Dynes formula the Authors have qualitatively re- 
constructed the experimental value of Tc- However, in 
the examined case the physical value of the Coulomb 
pseudopotential has been strongly lowered (^J. ^ 0.15). 

In the last step, we boldly underline that in the litera- 
ture exist other calculations of A than has been presented 



in the paper [T^]. In particular, Lei et al. have suggested 
a very large value of the electron-phonon coupling con- 
stant (A = 3.75 for p = 155 GPa and the sc structure) 
37| . On the other hand, Aftabuzzaman and Islam have 
predicted A = 0.903 for p = 161 GPa and the Pnma 
structure [l^- The last result is similar to the result 
obtained by Yin et al. [iTj . 



V. THE CHARACTERISTICS OF THE 
SOLUTIONS ON THE IMAGINARY AXIS 

The form of the order parameter on the imaginary axis 
for selected values of the temperature has been presented 
in Figure m (A). It has been shown that together with 
the increase of ujm the absolute values of A™ are de- 
creasing and are subjected to saturation. It should be 
underlined that taking the negative values by the order 
parameter function is connected with the non-zero value 
of the Coulomb pseudopotential. 

When analyzing the temperature's dependence of the 
order parameter, we found that absolute values of the 
function decrease together with the temperature's 
growth. The above result means that together with the 
growth of the temperature, the less number of Matsubara 
frequencies give significant contribution to the Eliashberg 
equations. 

The full dependence of the maximum value of the 
order parameter (A,„=i) on the temperature has been 
plotted in Figure |4] (B). We can observe that the val- 
ues 2Am=i {T) with a good approximation reproduce the 
temperature dependence of the energy gap at the Fermi 
level. 
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FIG. 4: (A) The dependence of the order parameter on Um for 
selected values of the temperature. (B) The maximum value 
of the order parameter as a function of the temperature. 

In Figure [5] (A) we have presented the form of the 
wave function renormalization factor on the imaginary 
axis. Similarly as for the order parameter, the increase 
of a;„i causes the decrease of the successive values of the 
function Zm- In the case of the high values of ujrm the 
function Zm is subjected to the saturation and takes the 
value equal to one. 
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Further, Figure [5] (B) presents the fuh dependence of 
the maximum value of the wave function renormahza- 
tion factor on the temperature. It can be noted that 
the presented function with a good approximation deter- 
mines the temperature dependence of the electron effec- 
tive mass. Moreover, from the obtained results, we can 
conclude that the electron effective mass takes a high 
value in the entire area of the existence of the supercon- 
ducting state. 
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FIG. 5: (A) The dependence of the wave function renormal- 
ization factor on ujm for selected values of the temperature. 
(B) The maximum value of the wave function renormalization 
factor as a function of the temperature. 



VI. THE PHYSICAL VALUE OF THE ORDER 
PARAMETER 

In order to determine the physical value of the order 
parameter for the chosen temperature, the solutions of 
the Eliashberg equations on the imaginary axis (?Wn) 
should be analytically continued on the real axis [ui). In 
the presented paper we have used the method introduced 
by Beach et al. [38| . The form of the order parameter on 
the real axis is being reproduced by using the function: 



PAl + •■• +PAr'^'' 



qAl + qA2LU 



QAr^ 



(10) 



where pAj and gAj denote the number coefficients, and 
r = 550. 

The dependence of the real and imaginary part of the 
order parameter on the frequency for selected values of 
the temperature has been presented in Figure |6l Addi- 
tionally, the rcscalcd Eliashberg function has been spec- 
ified. On the basis of the presented results, one can ob- 
serve that in the range of the low frequencies (from to 
about 20 meV), only the real part of the order parame- 
ter takes the non-zero values. From the physical point of 
view, the obtained result defines the lack of the damp- 
ing effects. For the higher values of the frequency (from 
about 20 mcV to about 40 meV), the real part of the 
order parameter takes relatively high values, which are 



clearly induced by the characteristic peaks in the Eliash- 
berg function. Furthermore, we can notice that in the 
discussed range of the energy, the imaginary part of the 
order parameter becomes non-zero and strongly increases 
together with the increase of frequency. For the higher 
frequencies (above 40 meV), the real part of the order 
parameter begins to vanish. This fact is related to the 
extinction of the Eliashberg function itself. 
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FIG. 6: The real and imaginary part of the order parameter 
on the real axis for selected values of the temperature. The 
rescaled Eliashberg function has been also presented. 

The physical value of the order parameter for the cho- 
sen temperature should be determined on the basis of the 
expression [l^, ps} : 



A (T) = Re [A (cj = A (T) , T)] 



(11) 



In the case of the superconductors the most interesting 
is the value of the order parameter for the temperature 
of zero Kelvin (A (0) ~ A (To))- On the basis of the sim- 
ple calculations we have made the following estimation: 
A (0) = 4.32 meV. 

Let us mention that familiarity with the value of A (0) 
and Tc allows to calculate the dimcnsionlcss ratio: Ri = 
2A (0) /ksTc- In the case of calcium we have obtained 
the following: 



Ri = 4.01. 



(12) 



The above result indicates that Ri considerably exceeds 
the value predicted by the BCS theory: [Ri] 

fH. 



BCS 



= 3.53 



VII. THE ELECTRON EFFECTIVE MASS 

The influence of the electron-phonon interaction on 
the electron effective mass (m*) can be determined on 



6 



the basis of the expression: m* = Re \Z (0)] me , where 
the symbol Z (0) denotes the value of the wave function 
renormalization factor on the real axis and rrif. is the bare 
electron mass. The form of the wave function renormal- 
ization factor on the real axis has been calculated with 
the use of the analytical continuation method: 



qzi + 



(13) 



where pzj and qzj arc the number coefficients, and r = 
550. 

In Figure [7] we have presented the shape of the func- 
tion Re[Z (w)] and lva[Z (w)] for the critical temperature. 
Similarly to the situation which took place in the case 
of the order parameter, for the low frequencies the non- 
zero is only the real part of the wave function renormal- 
ization factor. In the energy range around 20 meV we 
can observe characteristic but not so strong amplification 
of Re[Z'(a;)], which is clearly correlated with the peaks 
of the Eliashberg function. Additionally, the function 
Im[Z [ui)] is non-zero. In the range of the high frequen- 
cies, Re[Z (w)] decreases together with the increase of w. 
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FIG. 7: The real and imaginary part of the wave function 
renormalization factor on the real axis. Additionally, the 
rescaled Eliashberg function has been outlined. The inset 
represents the dependence of m^/me on the temperature. 

Next, the dependence of the ratio m*/me on the tem- 
perature has been determined. The results have been 
presented in the inset in Figure [T] We have found that 
the electron effective mass is large in the entire range, in 
which the superconducting state exists, and reaches its 
maximum equal to 2.32 for T = Tc- 

We can notice that for T = Tc the value of the ratio 
m*/me can be calculated with a great accuracy by us- 
ing the simple formula: m*/me ~ 1 -|- A = 2.27. The 
consistency between the exact numerical result and the 
analytical approach is the measure of the presented anal- 
ysis. 



From the physical point of view, the result presented 
above is particularly important, since it can be verified 
in a simple way if the measurement of the Sommerfeld 
coefficient is to be made in the future. 



VIII. THE THERMODYNAMIC CRITICAL 
FIELD AND THE SPECIFIC HEAT 

The thermodynamic critical field {Hq) and the differ- 
ence between the specific heat in the superconducting 
and normal state (AC = — C^) can be calculated on 
the basis of the free energy difference (AF = — F^): 
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(14) 
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The dependence of the free energy difference on the 
temperature has been presented in Figure [8] We can 
see that in the whole range of the existence of the su- 
perconducting phase, the value of the ratio AF/ p (0) is 
negative. From the physical point of view, it means that 
the superconducting state is thermodynamically stable. 

The thermodynamic critical field should be determined 
on the basis of the expression: 



He 



= V-Stt [AF/piO)]. 



(15) 



The influence of the temperature on the value of the ratio 
He I \/ P (0) has been presented in the Figure [H 
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FIG. 8: The ratios AF/p (0) and Hc/\/pW) as a function of 
the temperature. 

The difference of the specific heat has been determined 
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on the basis of the formula: 

AC _ l d^AF/p{0)] 
kBp{0) /3 dikeTf ■ ^ ' 

Additionally, the values of the specific heat in the normal 
state have been also determined: 

where 7 = Itt^ (1 + A). 

In Figure ^ we have plotted the dependence of the 
specific heat in the superconducting state and the normal 
state on the temperature. The characteristic "jump" , 
which appears at the critical temperature, can be easily 
noticed. 
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Taking into account the results above, we can state 
that the values of the considered ratios significantly di- 
verge from the values predicted by the classical BCS the- 
ory. In particular: [i?2]Bcs ~ ^■'^^ and [i?3]Bcg ^ 0.168. 

IX. SUMMARY 

In the paper, we have determined all relevant ther- 
modynamic parameters of the superconducting state in 
calcium under the pressure at 161 GPa. Wc have con- 
ducted all numerical calculations in the framework of the 
Eliashbcrg formalism, where the electron-phonon spec- 
tral function a^F (£7) has been taken form the paper 
On the basis of the exact numerical results, we can state 
that the depairing electron correlations in calcium are rel- 
atively strong (/i^ = 0.24). In the next step, the values 
of the parameters in the Allcn-Dyncs formula have been 
calculated. It has been shown that the critical tempera- 
ture is properly determined by the modified Allen-Dynes 
expression. 

Furthermore, we have proven that the thermodynamic 
properties of the superconducting state significantly di- 
verge from the predictions based on the simple BCS the- 
ory. In particular, the following values of the thermody- 
namic ratios have been obtained: i?i = 4.01, i?2 = 2.17, 
and i?3 = 0.158. 

In the last step, we have shown that the electron ef- 
fective mass is large in the entire area of the existence 
of the superconducting state, and [m*^^^^ = 2.32me at 
T = Tc. 
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